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Executive  Summary 


The  DMSP  Sensor  Fusion  Auroral  E-layer  Algorithm  builds  on  the  methods  used  in  the 
SSUSI  single  sensor  algorithm,  but  makes  use  of  information  from  the  SSJ/5  sensor  to  improve 
the  accuracy  of  the  deduced  E-layer  parameters.  The  Algorithm  operates  in  three  stages:  SSJ/5 
data  analysis,  Sensor  Fusion  analysis  of  subsatellite  SSUSI  pixels,  and  SSUSI  image  analysis  of 
off-track  pixels.  The  SSJ/5  data  analysis  characterizes  the  precipitating  electrons  and  protons 
along  the  satellite  track.  The  Sensor  fiision  analysis  combines  the  SSJ/5  data  analysis  with  the 
optical  intensities  observed  by  SSUSI  along  the  orbital  track  to  deduce  information  on  the  neutral 
composition,  in  particular  the  ratio  of  column  densities  of  atomic  and  molecular  oxygen  to  the 
column  density  of  molecular  nitrogen.  This  information  is  extrapolated  to  the  rest  of  the  image 
where  it  is  used  along  with  the  SSUSI  image  data  to  deduce  the  precipitation  characteristics  over 
the  whole  image.  This  information  is  then  used  to  calculate  the  peak  E-layer  electron  density  and 
the  height  of  the  E-layer  over  the  auroral  portions  of  the  image. 


V 


1.  Introduction 


The  ability  to  monitor  the  auroral  region  has  importance  for  satellite  operations, 
communications,  spacecraft  tracking,  and  optical  surveillance.  The  auroral  region  is 
characterized  by  the  presence  of  strong  fluxes  of  both  electrons  and  ions  with  energies  of  several 
keV.  These  particles  can  collect  on  the  surfaces  of  satellites  and  often  have  consequences  for  the 
operation  of  those  satellites.  These  particles  also  produce  auroral  optical  emissions  that  can 
interfere  with  ground  based  and  satellite  based  optical  surveillance  systems.  They  also  produce 
enhanced  ionization  in  the  E-layer  and  lower  F-layer  that  may  interfere  with  communications 
systems  or  with  radar  operations. 

Auroral  electrons  and  ions  can  be  detected  directly  only  by  flying  appropriate 
measurement  devices  through  the  precipitation  region  on  satellites.  To  map  the  region  even  at  a 
longitude  resolution  of  15°  would  require  24  satellites.  Since  it  is  not  practical  to  operate  more 
than  two  or  three  satellites  at  a  time,  mapping  the  precipitating  particles  must  be  done  using 
remote  sensing  of  the  effects  of  the  precipitating  particles  on  the  atmosphere. 

Auroral  optical  emissions  are  produced  by  the  impact  of  precipitating  electrons  and  ions 
with  energies  of  a  few  keV.  The  collisions  of  the  electrons  and  ions  with  the  constituents  of  the 
neutral  atmosphere  produce  atoms  and  molecules  in  excited  states.  These  excited  atoms  and 
molecules  dispose  of  their  excess  energy  by  radiating  at  wavelengths  ranging  from  infrared  to  X- 
rays. 


Because  the  atmosphere  is  largely  transparent  to  visible  and  infrared  emissions,  space 
based  observations  of  auroral  emissions  at  these  wavelengths  are  plagued  by  albedo  problems. 
Therefore,  the  ultraviolet  (UV)  wavelength  range  is  preferable.  The  relationship  between  particle 
precipitation  characteristics  and  UV  emissions  has  been  worked  out  using  electron  and  proton 
transport  models  [Strickland  et  al,  1983,  1984,  1993;  Daniell  et  ai,  1985;  Daniell  and 
Strickland,  1986].  These  papers  also  examined  the  influence  of  neutral  atmosphere  composition 
on  the  emissions  and  established  the  basic  principles  of  deducing  particle  precipitation 
characteristics  from  UV  measurements. 

The  single  sensor  SSUSI  Auroral  E-region  algorithm  (SSUSIAE)  has  been  described  in 
the  interim  Scientific  Report  No.  1  [Daniell  et  ai,  1994]  for  this  contract.  It  characterizes 
electron  and  proton  precipitation  in  the  diffuse  aurora  using  selected  ultraviolet  (UV)  emissions 
produced  in  the  E-region  of  the  auroral  ionosphere.  SSUSIAE  also  calculates  (the  peak 
electron  density  of  the  E-region)  and  h,„E  (the  height  of  the  E-region)  based  on  the  deduced 
particle  precipitation  characteristics.  It  is  limited  by  the  fact  that  the  effects  of  neutral 
composition  cannot  be  completely  eliminated,  introducing  a  large  uncertainty  in  the  deduced 
precipitation  characteristics. 

During  the  second  half  of  this  contract,  we  have  generalized  the  approach  used  in  the 
single  sensor  algorithm  to  incorporate  data  from  the  SSJ/5  instrument,  which  measures 
precipitating  electrons  and  ions  at  the  satellite.  The  incorporation  of  this  data  into  the  algorithm 
has  several  advantages.  First,  it  permits  better  characterization  of  the  precipitating  particles. 


1 


Second,  it  permits  the  determination  of  neutral  composition  information.  While  this  information 
is  strictly  available  along  the  satellite  track,  it  can  be  extrapolated  to  offtrack  pixels  with  varying 
levels  of  confidence. 

This  report  provides  a  ftinctional  description  of  the  algorithm  that  contains  sufficient 
detail  to  be  the  basis  of  a  Language  Independent  Description  from  which  a  full  implementation 
could  be  produced.  However,  it  should  be  noted  that  only  parts  of  the  algorithm  have  actually 
been  implemented,  and  others  have  not  been  defined  in  complete  detail.  This  description  should, 
therefore,  be  considered  provisional. 


2.  Top  Level  Description 

The  algorithm  consists  of  four  major  modules:  (1)  SSJ/5  data  analysis,  (2)  subsatellite 
sensor  fusion  analysis,  (3)  off-track  SSUSI  image  analysis,  and  (4)  the  Auroral  E-layer  analysis. 
The  first  three  modules  incorporate  new  analysis  techniques  that  have  not  been  included  in 
previous  sensor  algorithms  for  these  two  instruments.  The  basic  methods  used  in  each  module 
are  described  in  the  following  subsections.  The  algorithm  assumes  that  the  dayglow  and 
geocorona  contributions  have  already  been  removed  from  the  SSUSI  images. 


2. 1  SSJ/5  Data  Analysis 


The  first  module,  SSJ/5  data  analysis,  examines  electron  and  ion  spectra  to  determine  the 
precipitation  characteristics.  For  ions,  this  process  is  straightforward,  since  the  ion  spectra  are 
usually  well  characterized  by  a  Maxwellian  (thermal)  distribution  with  a  characteristic  energy  of 
a  few  keV.  The  general  form  of  a  Maxwellian  distribution  is 
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where  /  is  the  number  of  particles  with  energy  in  the  range  E  to  E+dE  and  velocity  vector  in  the 
solid  angle  represented  by  the  angular  ranges  0  to  0  +  c/0  and  cp  to  cp  -i-  c/(p .  E^  is  the 
characteristic  energy  (or  temperature)  of  the  distribution. 

For  a  Maxwellian  distribution,  the  differential  flux  of  particles  in  a  particular  direction, 
denoted  by  the  unit  velocity  vector  v  is 
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In  these  expressions,  Qm  is  the  downward  energy  flux,  obtained  by  integrating  z  ct)^  over  the 


,n 


downward  hemisphere  (y  <  0  <  rt ,  0  <  tp  <  27t ). 


The  shape  of  a  Maxwellian  distribution  is  completely  determined  by  a  single  parameter, 
Em-  The  peak  of  the  number  flux  (the  most  probably  energy)  is  Em,  the  mean  energy  E^  2Em, 
while  the  full  width  at  half  maximum  (FWHM)  is  ~2.9\Em.  The  peak  of  the  energy  flux  is  2Em. 
The  simplest  way  to  analyze  a  Maxwellian  spectrum  is  to  divide  the  differential  flux  by  E  and  to 
note  that 
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so  that  Em  is  the  inverse  slope  of  the  resulting  straight  line  and  Qm  may  be  obtained  from  the 
intercept.  This  is  best  done  using  the  highest  energy  channels  for  which  the  measured  flux  is 
well  above  the  instrument  sensitivity.  This  reduces  the  chance  of  error  due  to  the  presence  of  a 
background  of  low  energy  particles. 

For  electrons,  however,  the  spectrum  is  often  a  modifed  Maxwellian.  The  usual 
modification  is  a  shift  to  higher  energies  produced  by  parallel  electric  fields  [Hoffman,  1993], 
although  other  modifications  are  possible.  The  electron  data  analysis  distinguishes  three  cases 
(1)  the  diffuse  aurora,  characterized  by  a  thermal  (non-shifted)  Maxwellian,  (2)  an  inverted  V 
event  characterized  by  a  shifted  Maxwellian,  and  (3)  a  discrete  arc  characterized  by  nearly 
monoenergetic,  field  aligned  beams  that  often  deviate  substantially  from  a  simple  Maxwellian 
shape. 


The  basic  form  of  a  shifted  Maxwellian  is 
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As  for  the  unmodified  Maxwellian,  the  characteristic  energy  of  the  parent  Maxwellian  may  be 

(  1 

obtained  from  the  inverse  slope  of  the  straight  line  In 


■C 


However,  the  peak  energy  is  no 


longer  necessarily  determined  by  the  characteristic  energy  Em-  The  peak  of  the  spectrum  is  more 
pronounced  in  the  energy  flux,  where  the  peak  occurs  at  max(£’^,£:5).  Therefore, 

the  electron  analysis  must  determine  whether  the  primary  spectrum  is  shifted  or  not.  Unshifted 
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Maxwellians  are  assumed  to  be  isotropic  over  the  downward  hemisphere.  Shifted  Maxwellians 
are  usually  isotropic  over  most  of  the  downward  hemisphere. 

Even  after  the  primary  spectrum  has  been  characterized,  the  algorithm  must  check  for  a 
secondary  peak  somewhat  less  energetic  than  the  primary  peak.  This  peak  is  produced  by 
degraded,  backscattered  primaries  that  are  reflected  by  the  accelerating  potential  [Hoffman, 
1993].  These  secondary  peaks  are  usually  nearly  field  aligned  with  pitch  angles  less  than  30°. 
This  is  also  true  of  at  least  some  primary  peaks  in  intense  discrete  arcs  [McFadden  et  al,  1990]. 

The  electron  data  analysis  assigns  each  electron  spectrum  (obtained  once  per  second)  to 
one  of  the  three  classes  described  above.  Since  the  SSJ/5  instrument  does  not  provide  pitchangle 
information,  a  standard  pitch  angle  distribution  (different  for  each  class)  will  be  applied  to  each 
class.  The  downward  electron  energy  flux  (g^J  data  and  the  characteristic  electron  energy 
will  be  used  to  assign  a  crude  threshold  for  distinguishing  diffuse  aurora  from  inverted  V  and 
discrete  arcs.  These  will  be  used  in  the  SSUSI  image  analysis  to  assign  an  electron  spectrum 
classification  to  each  off  track  pixel. 


2.2  The  Sensor  Fusion  Analysis  for  subsatellite  pixels 

The  algorithm  must  trace  each  spectrum  along  the  spacecraft  field  line  down  to  an 
altitude  of  1 10  km  and  identify  the  pixel  containing  the  foot  of  the  field  line.  In  most  cases,  more 
than  one  spectrum  will  be  assigned  to  each  pixel.  The  electron  and  ion  spectral  characteristics 
(as  opposed  to  the  electron  and  ion  spectra  themselves)  will  be  averaged  over  the  pixel,  taking 
proper  account  not  to  blindly  average  different  classes  of  spectra.  The  differential  energy  flux 
will  be  used  to  assign  a  dominant  electron  spectrum  class  to  each  pixel. 

The  electron  and  ion  characteristics  deduced  from  the  SSJ/5  data  analysis  are  used  along 
with  121.6  nm,  135.6  nm,  LBHl,  and  LBH2  yield  curves  to  obtain  the  expected  values  of  the 
intensities  of  those  emission  features  using  a  nominal  neutral  atmosphere.  The  expected  values 
are  compared  to  the  observed  intensities  to  obtain  values  for  the  deviation  of  O/N2  and  O2/N2 
ratios  from  their  nominal  values. 

A  boundary  detection  algorithm  detects  the  equatorward  and  poleward  edges  of  the 
aurora.  Within  the  boundaries,  a  generalized  coordinate,  x,  is  defined.  If  the  equatorward 
latitude  is  and  the  poleward  latitude  is  then 


X  = 


(2.6) 


The  O/N2  and  O2/N2  ratios  are  treated  as  functions  ofjc.  Generally,  there  will  be  two  crossings  of 
the  auroral  oval  at  two  separate  locat  times  by  the  same  spacecraft.  Let  fo{x,f)  be  the  O/N2 
ratio  determined  at  the  first  local  time  and  fo{x,Tf)  be  the  ratio  determined  at  the  second  locat 
time.  Then  the  function  at  arbitrary  local  time  T  is 
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where 


The  only  requirement  being  that  A  similar  procedure  is  applied  to  the  O2/N2  ratio,  and 

the  ion  characteristic  energy.  (All  ions  are  assumed  to  be  protons.) 


2.3  The  SSUSI  image  data  analysis. 

For  each  pixel,  the  O/N2  and  O2/N2  ratios  are  determined  by  evaluating  the  appropriate 
function  of  the  form  given  by  Equation  (2.7).  The  analysis  now  proceeds  in  two  stages.  In  the 
first  stage,  the  analysis  assumes  that  both  protons  and  electrons  are  present.  Then  there  are  four 
parameters  {E^^,  E^p,  Q^p)  to  be  determined  from  four  data  elements  (/1216,  /i356,  /lbhi, 

/lbh2)- 


In  the  first  stage  analysis,  the  electron  spectrum  is  assumed  to  be  Maxwellian.  The  four 
parameters  and  the  four  data  elements  are  formally  related  by  the  four  (nonlinear)  equations 


^1216  “  2oe^216(^0’-^2»^0e)'*'2op^2l6(^’'^2  >^0p)  (2.8a) 

7lbhi  (2.8c) 

/lB„3  =e..jiBl,3(«oA,;£.,)  +  a,liB'i2(«0.«0,;-E.p)  <2.8d) 


where  the  functions  etc.  represent  the  yield  of  the  indicated  emission  from  a  downward 
energy  flux  of  1  erg  cm"^  s  *  carried  by  the  indicated  particles.  (Of  course,  1^2  le  vanishes 
identically.)  Rq  (Rq^)  is  the  scale  factor  for  the  atomic  oxygen  density  (molecular  oxygen 
density)  from  a  model  (e.g.,  MSIS-90)  required  to  make  the  model  give  an  O/N2  (O2/N2)  ratio 
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that  conforms  to  the  estimate  extrapolated  from  the  sensor  fusion  analysis.  If  and  are 
known,  then  we  have  four  (nonlinear)  equations  in  four  unknowns. 


The  Second  Stage  of  the  analysis  depends  on  the  results  of  the  first  stage  analysis. 


Case  1 :  If  0.05  <  <  20.0  then  both  protons  and  electrons  are  present.  Consider  two  cases: 

Qop 

Case  1  a:  If  then  the  precipitation  is  probably  due  to  the  diffuse  aurora,  so 

Accept  Stage  1  analysis  without  modification 
Case  lb:  If  >  Q„,re,hoid  the  precipitation  is  probably  due  to  an  inverted  V  event  or 
a  discrete  arc,  so  Redo  Stage  1  analysis  using  monoenergetic  electron  yield 
functions  rather  than  Maxwellian  electron  yield  functions. 

Case  2:  If  -^^<0.05,  then  electrons  are  effectively  absent.  Therefore,  perform  Pure  Proton 

Qop 

analysis:  4  equations,  4  unknowns: 

^216 
-^1356 
^LBHI 
^LBH2 


(2.9a) 

(2.9b) 

(2.9c) 

1  “  0Op^BH2(^’^O2  ’  ^Op) 

(2.9d) 

where  the  four  unknowns  are  R^y,  Rq  ,  ,  and  Q 


cap- 


Qoe 


Case  3:  If  — ^>20.0,  then  protons  are  effectively  absent.  Therefore,  perform  Pure  Electron 


a 


•Op 


analysis:  3  equations,  4  unknowns: 


'1356 


'LBHI 


'  LBH2 


=  Qoe  ^356  ( ■^  » -^2  ’  ) 

~  Goe  ^BH2  (  ^  ’  ^Oe  ) 


(2.10b) 

(2.10c) 

(2.10d) 


One  of  the  four  unknowns,  R^,  R^^,  Eq^,  and  cannot  be  determined  from  the  UV 

intensities  alone.  Since  O2  is  likely  to  be  less  variable  than  O,  the  analysis  assumes  that 
Rq^  is  given  by  Equation  (2.7).  The  remaining  three  unknowns  are  now  determined  by 

Equations  (2.10). 


2.4  The  Auroral  E~layer  Characteristics 

The  E-layer  characteristics,  f„E  and  h^E,  are  determined  in  much  the  same  way  as  in  the 
single  sensor  algorithm,  using  an  altitude  dependent  effective  recombination  coefficient. 
However,  in  the  sensor  fusion  algorithm,  the  effective  recombination  coefficient  will  be  a 
function  of  R^y  and  R^  . 
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3.  SSJ/5  Data  Analysis  Module 

This  module  ingests  SSJ/5  data,  analyses  each  ion  and  electron  energy  spectrum  to 
determine  the  precipitation  characteristics,  defined  as 

The  characteristic  (most  probable  or  peak)  energy  in  keV  of  the  ions  (assumed  to 
be  protons). 

0,^  The  downward  energy  flux  of  the  ions  (erg  cm'^  s') 

E^^  The  characteristic  (most  probable  or  peak)  energy  in  keV  of  the  main  electron 
spectrum. 

The  downward  energy  flux  of  the  main  electron  spectrum  (erg  cm'^  s  ') 

ST^  The  spectral  type  of  the  main  electron  spectrum  (Maxwellian,  shifted  Maxwellian, 
Monoenergetic). 

SP  A  flag  indicating  the  presence  or  absence  of  a  secondary  peak  in  the  electron 
spectrum.  If  present,  it  is  assumed  to  consist  of  degraded  primaries,  backscattered 
by  the  atmosphere  and  reflected  by  the  auroral  potential  region. 

£,^  The  characteristic  (most  probable  or  peak)  energy  in  keV  of  the  secondary 

electron  spectrum.  ^  ^ 

The  downward  energy  flux  of  the  secondary  electron  spectrum  (erg  cm'^  s'). 

The  ions  are  assumed  to  always  have  a  Maxwellian  distribution.  For  the  electrons,  note  that  the 
“secondary  electron  spectrum”  is  different  than  the  “spectrum  of  the  secondary  electrons.” 


3.  /  SSJ/5  Data  Ingestion  and  processing 

The  SSJ/5  data  comes  in  the  form  of  differential  number  fluxes  of  ions  and  electrons  in 
20  energy  channels  ranging  from  30  eV  to  30  keV.  The  channels  are  organized  into  two  groups 
of  10.  Channels  1  through  10  represent  energies  from  30  keV  to  1  keV.  Channels  1 1  through  20 
represent  energies  from  1  keV  to  30  eV.  Channels  10  and  11  both  measure  1  keV  particles. 
Since  the  instrument  steps  through  the  two  groups  of  channels  in  parallel,  Channel  10  is  counted 
at  the  end  of  the  sweep  while  Channel  11  is  counted  at  the  beginning  of  the  sweep. 
Consequently,  the  two  measurements  of  1  keV  particles  are  often  slightly  different.  A 
comparison  of  these  two  channels  gives  an  indication  of  the  amount  of  temporal  aliasing  present 
in  the  data. 

A  preprocessor  for  SSJ/4  data  was  provided  for  use  with  the  PRISM  algorithm.  The 
PRISM  preprocessor  provides  only  mean  energy  (keV)  and  energy  flux  (erg  cm  ^  s  ')  The  AEL 
Algorithm  requires  more  extensive  processing  of  the  SSJ/5  data  because  it  needs  to  distinguish 
between  the  diffuse  aurora  (characterized  by  broad,  approximately  Maxwellian,  spectra)  and 
discrete  arcs  (characterized  by  narrow  spectra,  often  shifted  Maxwellians).  Therefore,  the  AEL 
Algorithm  assumes  the  existence  of  a  preprocessor  for  the  SSJ/5  data  whose  output  is  the 
differential  flux  in  each  energy  channel. 
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MODULE  1;  SSJ/5  data  ingestion 


PURPOSE:  To  ingest  preprocessed  SSJ/5  data,  correct  for  temporal  aliasing,  and  reverse  the 
energy  ordering  of  the  data. 

INPUT: 

8**^*  =  e{‘^*,8^‘^^...E^o  j,  the  vector  of  electron  energy  channels  (eV) 

g(p)  _  the  vector  of  ion  energy  channels  (eV) 

f  (f)  ^  the  vector  of  electron  differential  particle  fluxes  (cm*^  s  '  eV'  sr' 

') 

f(p)  ^  the  vector  of  ion  differential  particle  fluxes  (cm'^  s'  eV'  sr') 

s**^'  =  \-.-4o  ]» vector  of  electron  channel  sensitivities,  i.e.,  the  differential 

electron  flux  corresponding  to  1  count 

si/')  _  the  vector  of  ion  channel  sensitivities,  i.e.,  the  differential  ion  flux 

corresponding  to  1  count 

NOTE:  For  both  electrons  and  ions  e,  >82  >...e,o  and  e,,  >8,2  >...820.  £,0  and  8,,  are  both 
approximately  1000  eV,  but  no  special  ordering  of  these  two  energies  is  assumed. 

OUTPUT: 

the  vector  of  electron  channel  energies  with 

_  (|)|0  (jjW  the  corresponding  vector  of  electron  channel  fluxes  (cm’^  s''  eV' 

sr''). 

E*/')  =  j,  the  vector  of  ion  channel  energies  with  E^^j  >  E^^* 

(j,!/-)  _  (j,|^),(|)U)  ,j,W  the  corresponding  vector  of  ion  channel  fluxes  (cm'^  s  '  eV'  sr  '). 

the  vector  of  electron  channel  sensitivities,  i.e.,  the  differential 
electron  flux  corresponding  to  1  count 

S(/')  =  the  vector  of  ion  channel  sensitivities,  i.e.,  the  differential  ion 

flux  corresponding  to  1  count 

NOTE:  For  both  electrons  and  ions,  the  ordering  of  the  energy  channels  has  been  reversed 
relative  to  the  input  convention,  and  one  of  the  redundant  I  keV  channels  has  been 
removed. 

Step  1 :  Interpolate  channels  10  and  1 1  to  precisely  1  keV. 

Definition:  is  the  1  keV  electron  flux  inferred  from  the  channel  1 0  flux, 
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Definition:  is  the  1  keV  electron  flux  inferred  from  the  channel  1 1  flux, 

Definition:  is  the  1  keV  proton  flux  inferred  from  the  channel  10  flux, 

Definition:  F.\^^  is  the  1  keV  proton  flux  inferred  from  the  channel  1 1  flux, 


^.0  =  exp  \nflo’  + 


=  exp  ln/,P  + 


1000-eS? 

pW-pW 
^9  ^10 

1000 -el"’ 


'1  ^  pW_pW  ' 

t>12  t>|j 

innn_o(^) 


(in/,'" -In/,; 
(in/,;" -In/,; 


F<'’>-exn  In  .rfp)  ■  ^,0  r(p)_j^  Ap) 

^10  —  exp  in/,g  +  ^^iny,  inyio 


10  -^(P)_„(P) 

L  '’9  *^I0 

cxnln/-'".‘"””~^''* 
/<,,  -exp  my,,  +  (p)_  (p) 

fc|2  ^11 


(in/,'/' -In/,; 


Step  2:  Determine  the  change  in  1  keV  flux  over  the  duration  of  the  sweep: 

A/rlp) 

Step  3:  Refer  each  channel’s  flux  back  to  the  start  of  the  sweep: 


For/=  1,2, ...9  do 
y  =  20  -  i 

AF 

^  j  g  p{e) 

L  ^  ^10 

Q  pip) 
^10 

Fj"’=e|"’ 

=  ej'’’ 

=5*"’ 

=  4"^ 

End  do 

X(-)  ^  fie) 

Yio  “ 
x(p)  ^  f(p) 

9 10  ~Jli 

r(^)  _ 

^10  ■"*^11 

p(p)  — 

^10 

*-'10  ”*^11 
c(/^) 

Oio  —All 
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11,12, 

...19( 

do 

y  =  20 

-  i 

<l>r  = 

1-i 

-10 

9 

•^10 

Ji+l 

i-i 

-10 

9 

^ipy 

p{p) 

Ap) 

Jn-\ 

Ef  = 

^!:l 

E^f  = 

p(/’) 

*^/+i 

5)^*  = : 

= , 

End  do 

END  MODULE  1. 


3.2  Electron  data  analysis 

Electron  and  ion  data  are  analyzed  separately  because  the  two  types  of  auroral  particles 
are  subject  to  different  physical  processes.  The  algorithm  recognizes  three  types  of  electron 
spectra:  Maxwellian,  Shifted  Maxwellian,  and  Monoenergetic.  (Technically,  shifted 

Maxwellians  are  a  kind  of  monoenergetic  spectrum.) 

Maxwellian  spectra  are  characteristic  of  the  diffuse  aurora,  which  is  produced  by  a 
population  of  (hot)  thermal  electrons.  Such  a  spectrum  is  broad  and  isotropic.  It  is  often 
accompanied  by  a  low  energy  component  consisting  of  secondary  electrons.  This  low  energy 
component  has  little  effect  on  the  auroral  E-layer/ 

Shifted  Maxwellian  spectra  a  characteristic  of  inverted-V  events  and  consist  of  thermal 
electrons  that  have  been  accelerated  (and  sometimes  heated)  by  passage  through  a  field-aligned 
potential  drop.  The  spectral  shape  is  still  Maxwellian,  but  the  peak  energy  is  much  higher  than 
the  characteristic  energy  of  the  Maxwellian.  Since  the  width  of  the  spectrum  is  no  much  smaller 
than  the  peak  energy,  the  spectrum  appears  to  be  nearly  monoenergetic.  The  spectrum  is  usually 
nearly  isotropic.  A  secondary  peak  in  the  spectrum  is  often  present.  This  secondary  peak 
consists  of  degraded  primaries  that  are  backscattered  up  the  magnetic  field  line  but  are  reflected 
by  the  auroral  potential  region.  The  secondary  peak  is  nearly  field  aligned. 

Monoenergetic  spectra  a  characteristic  of  discrete  arcs  and  other  auroral  features  in  which 
the  Maxwellian  shape  of  the  parent  population  of  electrons  is  sufficiently  distorted  so  that  a 
Maxwellian  based  analysis  is  not  possible.  A  secondary  peak  and  other  confounding  structure 
may  or  may  not  be  present. 
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MODULE  2:  SSJ/5  electron  data  analysis 

PURPOSE:  To  determine  whether  the  electron  spectrum  may  be  characterized  as  Maxwellian, 
shifted  Maxwellian,  or  monoenergetic,  and  to  determine  the  characteristics  of  the 
main  and  (if  present)  secondary  peaks. 

INPUT: 

the  vector  of  electron  channel  energies  (eV)  with 
=  ([)**'*, the  corresponding  vector  of  electron  channel  fluxes  (cm*^  s'*  eV* 
sr'*). 

S**"’  =  the  corresponding  channel  sensitivity,  i.e.,  the  flux  corresponding 

to  1  count. 

OUTPUT: 

57^,,  a  string  variable  containing  “Maxwellian”  if  the  electron  flux  appears  to  be  an 

unshifted  Maxwellian,  or  “Shifted  Maxwellian”  if  the  electron  flux  appears  to 
be  a  shifted  Maxwellian,  or  “Monoenergetic”  if  the  electron  flux  appears  to  be  a 
narrow  peak  of  uncertain  or  non-Maxwellian  shape, 
the  characteristic  (maximum  or  most  probable)  energy  (keV)  of  the  main  peak  of 
the  electron  spectrum. 

00^,  the  downward  energy  flux  (erg  cm*^  s'*)  of  the  main  peak  of  the  electron  spectrum, 
the  characteristic  energy  of  the  parent  Maxwellian  when  the  electron  spectrum  is 
identified  as  a  “Shifted  Maxwellian.”  For  “Maxwellian”  or  “Monoenergetic” 
spectra,  this  variable  is  exactly  0. 

a  logical  variable  with  the  value  TRUE  if  a  secondary  electron  peak  is  detected, 
and  the  value  FALSE  otherwise. 

E,^,  the  characteristic  energy  (keV)  of  the  secondary  peak  (if  any)  of  the  electron 
spectrum.  If  SECONDARY  is  FALSE,  this  variable  is  undefined, 
the  downward  energy  flux  (erg  cm'^  s  ')  of  the  secondary  peak  (if  any)  of  the 
electron  spectrum.  If  SECONDARY  is  FALSE,  this  variable  is  undefined, 
the  full  width  at  half  maximum  (FWHM)  of  the  secondary  peak  (if  any)  of  the 
electron  spectrum.  If  SECONDARY  is  FALSE,  this  variable  is  undefined. 

Begin  Procedure 

Step  1 .  Locate  the  absolute  maximum  in  the  energy  flux  spectrum: 

Definition:  is  a  vector  containing  the  energy  flux  spectrum  (eV  cm'^  s'*  sr  '). 

Definition:  is  the  subscript  or  index  (channel  number)  of  the  maximum  element  of 

(e) 

\|/^  ^ 

Definition:  \|/  is  the  value  of  the  maximum  element  of  (eV  cm'^  s''  sf'). 

Definition:  is  the  energy  (eV)  of  channel 
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Definition;  MAX  is  a  function  procedure  or  module  that  accepts  as  input  a  vector  of 
numerical  values  and  returns  the  largest  of  these  along  with  it  subscript. 


Step  2.  Determine  the  characteristics  of  the  primary  spectrum. 


Definition: 

Definition: 


Definition: 

Definition: 

Definition: 

Definition: 


is  the  phase  space  flux  density 

LLSQ  is  a  procedure  or  module  that  accepts  a  vector  ofjc  values  and  a  vector 
of  y  values  and  performs  a  linear  least  square  fit  of  the  function  ax  +  b  io 
those  values. 

is  the  number  of  points  from  the  electron  spectrum  to  be  used  in  LLSQ 
is  the  characteristic  energy  (eV)  of  a  Maxwellian  energy  spectrum. 

E^  is  the  energy  shift  (eV)  of  a  Shifted  Maxwellian  energy  spectrum. 

is  the  downward  energy  flux  (eV  cm'^  s' ‘)  of  a  Maxwellian  energy 
spectrum  that  is  isotropic  over  the  downward  hemisphere 


=  19 

j  =  /max  +  1 

While/„„^  =  19do 
If  (j)!;!  <  Then 

J  stop  J  1 

End  If 

y=y+i 

End  While 

If  ^ymax  +  1  Then 

The  spectrum  cannot  be  analyzed  in  terms  of  Maxwellians 
Else 

N  =  i  -  / 

fit  J  stop  J  max 

End  IF 


n 


(‘■)  _ 


f,"'  (1) 


7(‘')  ’  rie) 


TI9 

rW 

£-19 


If  Mji,  =  2  Then 

The  two  points  define  a  straight  line: 
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a  = 


*  Jm&\ 


(e) 

■/nnx+l 


E(‘) 

^y„»x+2 


■^7n»,+l 


Z)  = 


y™x+2'L/m5x+l 


^7'J,+I^Lx+2 


,(«) 


f(^)  _  f(^) 

7inax+2  ^max"^^ 


Else 


Use  LLSQ 


y  = 


X  = 


r(^) 

■^7max+' 


r(^)  F 

’  ymax  ’  Ato/? 


r(^) 


LLSQ(iV 
End  If 


f" 


,x,y. 


,a,i)) 


NOTE:  For  a  Maxwellian,  the  peak  of  the  energy  flux  occurs  at  twice  the  characteristic 
energy.  Thus,  we  expect  for  an  unshifted  Maxwellian.  If  this  condition 

holds,  then  we  classify  the  spectrum  as  diffuse  aurora.  If  E^  «  ,  then  we 

classify  the  spectrum  as  a  discrete  arc. 


5£’„,ax  = 

=  2E^ 

Emm 

A,  = 

- 

^  low 

ynm 

7.mx-l 

^high  - 

■ 

■  ^.w+1 

1 

. — . 

i 

If6£„„iOThen 

If|5E,,,J<A,,„Then 

The  data  are  consistent  with  a  Maxwellian  energy  spectrum 
ST^  =  "Maxwellian” 

Eoe  ~ 

Qm  ~  ^kE/^c 

a,=i.6xio-'’e„ 


Else 

The  data  are  consistent  with  a  shifted  Maxwellian  energy  spectrum 
ST^  =  "Shifted  Maxwellian" 


Es  ~  ^max 
Eo,  =  10-^E_ 

Ems  ~ 


Qm  ~  ^^{Em  E mEs  +  2  Es)E exp 


b- 


eo,  =  1.6xlO-'^xe^ 

End  If 
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Else  If  >  0  Then 
Then 

The  data  are  consistent  with  a  Maxwellian  energy  spectrum 
ST^  =  "  Maxwellian" 

=  2KEy 
a.,  =  l6xlO-'^Q^ 

Else 

The  data  are  not  consistent  with  either  a  Maxwellian  or  a  Shifted  Maxwellian 
energy  spectrum. 

Details  are  TBD. 

End  If 
End  If 


Step  3.  Model  the  primary  spectrum  according  to  the  characteristics  determined  in  Step  2. 

Definition:  is  the  model  primary  electron  spectrum  based  on  the  characteristics 

determined  in  Step  2. 

Definition:  MXW  is  a  function  procedure  that  calculates  the  differential  number  flux  of  a 
(possible  shifted)  Maxwellian  given  Q^,  and  Ef 


MXW(£,e„, = 


Cm 


exp 


2n{El  +  E^E,+iE^,)E^ 

0 

Definition:  is  the  maximum  value  of  the  model  flux: 

^  _|MXW(£„,e„,£„£j),  £„>£ 

|MXW(£„e„, £„£,),  £„<£ 
Definition:  Eq^^^  is  the  energy  corresponding  to 


'  M  J 


,  E  >  E. 


E  <  E^ 


4i'''=[MXW(£,''',e«.£„,£s),MXW(£i'',e„,£„.£s)....MXW(£;,'',e„,£„,£,) 


If  Ej  <  E^  Then 


a 


M 


2n{El*E„E,+\El] 


exp 


-1 


^Omax  ~  ^ M 


Else 


Cm 


2ji(£j, +£„£,+!£')£ 


M 


End  If 


Step  4.  Determine  the  characteristics  of  the  secondary  peak  (if  it  is  present). 


14 


Definition:  is  the  measured  differential  electron  spectrum  after  subtraction  of  the 

modeled  primary  spectrum  from  Step  3. 

Definition:  is  the  maximum  element  of 

Definition:  is  the  subscript  or  index  of  ^^3^ 

Definition:  £,  is  the  energy  of  channel 
Definition:  co  ^3^  is  the  differential  energy  flux  at  £■, 

£,  = 

^  max  —  ^l^max 

Iff0„3,  <10-'£omax^maxThen 

The  flux  in  the  secondary  peak  is  too  small  to  matter 

^\e  ~2^0e 

e,.=o 

Else 

The  flux  in  the  secondary  peak  should  not  be  neglected 
Find  FWHM  -  details  TBD. 

End  If 


END  MODULE  2. 


3.3  Proton  data  analysis 

Electron  and  ion  data  are  analyzed  separately  because  the  two  types  of  auroral  particles 
are  subject  to  different  physical  processes.  The  proton  algorithm  assumes  that  the  energy 
spectrum  is  approximately  Maxwellian. 

MODULE  3:  SSJ/5  proton  data  analysis 

PURPOSE:  To  determine  the  characteristics  of  the  proton  energy  spectrum. 

INPUT: 

E^p)  =  the  vector  of  proton  channel  energies  (eV)  with 

the  corresponding  vector  of  proton  channel  fluxes  (cm'^  s''  eV' 
sr''). 
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the  corresponding  channel  sensitivity,  i.e.,  the  flux 
corresponding  to  1  count. 

OUTPUT: 

the  characteristic  (maximum  or  most  probable)  energy  (keV)  of  the  main  peak  of 
the  proton  spectrum. 

00^,  the  downward  energy  flux  (erg  cm'^  s  ’)  of  the  main  peak  of  the  proton  spectrum. 


Begin  Procedure 


Step  1 .  Locate  the  absolute  maximum  in  the  energy  flux  spectrum: 


Definition: 

Definition: 

Definition: 

Definition: 

Definition: 


is  a  vector  containing  the  energy  flux  spectrum  (eV  cm'^  s"'  sf'). 
is  the  subscript  or  index  (channel  number)  of  the  maximum  element  of 

\|/  is  the  value  of  the  maximum  element  of  (eV  cm'^  s'  sf'). 

^max  is  the  energy  (eV)  of  channel 

MAX  is  a  function  procedure  or  module  that  accepts  as  input  a  vector  of 
numerical  values  and  returns  the  largest  of  these  along  with  it  subscript. 


v™„=MAX(v'‘'',y,„) 


Step  2.  Determine  the  characteristics  of  the  proton  spectrum. 


Definition:  is  the  phase  space  flux  density 

Definition:  LLSQ  is  a  procedure  or  module  that  accepts  a  vector  ofx  values  and  a  vector 
of  y  values  and  performs  a  linear  least  square  fit  of  the  function  ax  +  b  to 
those  values. 

Definition:  Njj,  is  the  number  of  points  from  the  proton  spectrum  to  be  used  in  LLSQ 

Definition:  is  the  characteristic  energy  (eV)  of  a  Maxwellian  energy  spectrum. 

Definition:  Qj^,  is  the  downward  energy  flux  (eV  cm'^  s  ')  of  a  Maxwellian  energy 
spectrum  that  is  isotropic  over  the  downward  hemisphere 

J  =  imax  +  1 

Whiley„„,  =  19do 
If  (f)*'’*  <  sj/’*  Then 

J  stop  J  1 
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End  If 
7 = y + 1 
End  While 

Ify.«„p^y..ax+lThen 

The  spectrum  cannot  be  analyzed  in  terms  of  Maxwellians 
Else 

^ fit  j stop  Jmax 

End  IF 


Ap)  _ 


♦I'’  V/' 


(P) 

19 


p{p)  ’  p(p)  ’  pip) 
^2  ^19 


If  =  2  Then 


The  two  points  define  a  straight  line: 

(P) 

'yimx+^ 


a  = 


•nTii!l.+2 


*  yniax ' 

?(P) 


rip) 


pypi  _ 

./ni8x  Ji\vu 


b  = 


./ max 


ip)  „(p) 


pip)  _  I7(P) 

./iiiax'*'^  -/max'*'^ 


Else 


Use  LLSQ 


y  = 


X  = 


Inri*/’*  ,lnr|*/*  ^2,...lnT|^^’ 

'./max+l  *  yniax  'Jstop  ^ 


iP)  pip)  pip)  1 

/max'*'^  ’  yniax"**^’  Jxtop  J 

LLSQ(Ar^,,x,y,a,^>) 


End  If 


E^=-- 

a 


NOTE:  For  a  Maxwellian,  the  peak  of  the  energy  flux  occurs  at  twice  the  characteristic 
energy.  Thus,  we  expect  Ej^  »jE^  for  an  unshifted  Maxwellian. 


A  =  -  F'.'' 

/max  yinax  ^ 

A  =  F*'"'  -  F*''* 

^Itigh  ^,/„„x+'  >imx 

If  5£.„  <  0  Then 

If|5£„„|<A,„.Then 

The  data  are  consistent  with  a  Maxwellian  energy  spectrum 

=2KEy 


rip) 


17 


Q,^  =  l6xlO-'^Q^ 

Else 

The  data  are  not  consistent  with  a  shifted  Maxwellian  energy  spectrum 
Return  Error  Condition  TBD. 

End  If 

Else  If  >  0  Then 
Then 

The  data  are  consistent  with  a  Maxwellian  energy  spectrum 
Eo,  =  10-^Em 

Qm  ~ 

Qo,  =  l6xlO-'^Q^ 

Else 

The  data  are  not  consistent  with  a  Maxwellian  energy  spectrum. 

Return  Error  Condition  TBD. 

End  If 
End  If 

END  MODULE  3. 


4.  SSUSI  On-Track  Data  Analysis  Module 

This  module  ingests  SSUSI  121.6  nm,  135.6  nm,  LBHl,  and  LBH2  images,  identifies  the 
pixels  that  lie  on  the  DMSP  orbital  track,  and  analyses  the  image  data  in  conjunction  with  the 
SSJ/4  data  (sensor  fusion).  The  SSUSI  data  is  assumed  to  consist  of  intensities  (Rayleighs) 
associated  with  each  pixel.  The  dayglow  and  geocoronal  backgrounds  are  assumed  to  have  been 
subtracted,  leaving  only  the  auroral  contribution.  Each  pixel  is  assumed  to  be  associated  with  a 
geographic  location  (latitude,  longitude),  geomagnetic  location,  and  solar  zenith  angle.  The 
notation  used  in  this  section  includes 

/,2|(, ;  The  1 2 1 .6  nm  (Lyman  alpha)  intensity  due  to  proton  precipitation.  (R) 

/i35(,:  The  135.6  nm  intensity  due  to  electron  and  ion  precipitation  (R) 

/lbhi  :  The  LBHl  (140-150  nm)  intensity  due  to  electron  and  ion  precipitation  (R) 

/lbh2‘  The  LBH2  (165-180  nm)  intensity  due  to  electron  and  ion  precipitation  (R) 

The  geographic  latitude  of  the  pixel  (deg) 

(pg :  The  geographic  longitude  (deg) 

:  The  geomagnetic  latitude  (deg) 

(p„, :  The  geomagnetic  longitude  (deg) 

X  :  The  solar  zenith  angle  (deg) 

T  ;  Universal  Time  (UT,  hours) 

T/.  Solar  Local  Time  (SLT,  hours) 
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:  Magnetic  Local  Time  (MLT,  hours) 


4. 1  SSUSl  Data  Ingestion  and  processing 

The  SSJ/5  data  comes  in  the  form  of  intensities  in  four  wavelength  channels.  The 
intensities  are  in  Rayleighs  (R). 

7,2,^  is  the  intensity  of  the  Lyman  alpha  line  of  atomic  hydrogen  at  122  nm.  The 
auroral  component  comes  exclusively  from  precipitating  protons  and  hydrogen 
atoms.  There  is  no  electron  component. 

7,35^  is  the  combined  intensity  of  the  atomic  oxygen  2/?“  -3s  transition  at 

136  nm  and  the  (3,0)  band  at  135  nm  of  the  Lyman-Birge-Hopfield  system  of 
molecular  Nitrogen.  Both  protons  and  electrons  may  contribute  to  the  observed 
intensity. 

7lbhi  is  the  intensity  of  all  of  the  LBH  bands  in  the  140-150  nm  channel,  along 
with  the  atomic  nitrogen  Ip^^D-  35  transition  at  149  nm.  Both  protons  and 
electrons  may  contribute  to  the  observed  intensity. 

7,0„2  is  the  intensity  of  all  of  the  LBH  bands  in  the  165-180  nm  channel,  along 
with  the  atomic  nitrogen  lp^^P-3s^P  transition  at  1 74  nm.  Both  protons  and 
electrons  may  contribute  to  the  observed  intensity. 

MODULE  1:  SSUSl  data  ingestion 

PURPOSE;  To  ingest  preprocessed  SSUSl  image  data.  The  precise  definition  of  this  module 
depends  on  the  precise  way  in  which  it  is  to  interact  with  the  SSUSl  single  sensor 
algorithm. 

INPUT; 

The  location  of  the  SSUSl  image  data. 

NOTE;  The  information  required  to  specify  the  location  of  the  SSUSl  data  depends  on  the 
details  of  the  SSUSl  single  sensor  algorithm. 


OUTPUT; 

n^,  the  number  of  pixels  in  the  cross-track  direction 
n, ,  the  number  of  pixels  in  the  along-track  direction 
Ii2i(,,  an  n,  by  n,  dimensional  array  containing  the  121.6  nm  intensity  7,2,6 
the  n^  X  n^.  pixels 
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I,35(,,  an  by  n^,  dimensional  array  containing  the  135.6  nm  intensity  for  each  of 
the  X  rt,  pixels 

Ilbhi,  an  by  tiy  dimensional  array  containing  the  LBHl  intensity  /lbh,  for  each  of  the 
«,  X  n^.  pixels 

Ilbh2?  an  by  dimensional  array  containing  the  LBH2  intensity  /lbh2  for  each  of  the 
X  pixels 

G,  an  by  n^.  dimensional  array  containing  the  geographic  coordinates  for  each  of  the 
X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

M,  an  by  n^.  dimensional  array  containing  the  geomagnetic  coordinates  for  each  of 
the  X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

X,  an  by  n^.  dimensional  array  containing  the  solar  zenith  angle  x  for  each  of  the 
n^xn^,  pixels. 

vj/,  an  by  dimensional  array  containing  the  nadir  angle  \j/  for  each  of  the  x  n^. 
pixels. 

NOTE;  Individual  elements  of  the  above  arrays  will  be  denoted  by  a  symbol  of  the  form 
A(m,n)  where  A  is  an  by  array,  m  represents  the  column  and  n  the  row  of  the 
particular  pixel.  See  Figure  1  for  an  illustration. 


Step  1 :  Access  image  data,  reformat  as  necessary,  store  in  output  arrays.  Details  TBD. 

END  MODULE  1. 
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Figure  1  SSUSI  Pixel  Identification  Scheme 

4.2  Identification  of  subsatellite  pixels 

The  electrons  and  ions  measured  by  the  SSJ/5  instrument  are  guided  by  field  lines  and 
deposit  their  energy  in  the  atmosphere  threaded  by  the  field  line  passing  through  the  satellite  at 
the  time  of  the  measurement.  Thus,  for  the  purposes  of  this  algorithm,  “subsatellite”  does  not 
indicate  the  pixel  directly  beneath  the  satellite,  but  the  pixel  through  which  the  satellite’s  field 
line  passes.  The  algorithm  assumes  that  the  emitting  layer  is  at  1 10  km,  so  it  traces  the  field  line 
down  from  the  satellite  altitude  to  that  altitude  to  determine  through  which  pixel  it  passes. 


MODULE  2:  Identification  of  subsatellite  pixels 

PURPOSE:  To  trace  field  lines  from  the  satellite  to  the  E-layer,  nominally  110  km,  and  to 
identify  subsatellite  pixels.  Note  that  these  will  not  necessarily  be  pixels  along 
the  central  column  of  the  image  because  the  algorithm  defines  subsatellite  as  lying 
on  the  same  field  line  as  the  satellite  but  at  an  altitude  of  1 10  km.  Since  there  will 
generally  be  more  than  one  spectra  corresponding  to  the  same  pixel,  the  module 
produces  an  “average”  set  of  spectral  characteristics  for  each  subsatellite  pixel. 

INPUT: 
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N^,  the  number  of  electron  and/or  ion  spectra  taken  during  the  same  time  interval  as 
the  image 

^  ^  1 2  ^  /V,  ]» the  vector  of  times  (UT)  corresponding  to  each  of  the  spectra 

X  =  I ,  X  2 , . . .  A,  ^  j,  the  vector  of  satellite  geographic  latitudes  corresponding  to  each  of 

the  spectra 

9  =  [9 1 5^2  ’  /V,  ]’  the  vector  of  satellite  geographic  longitudes  corresponding  to  each  of 

the  spectra 

h  =  ,^2  j,  the  vector  of  satellite  altitudes  corresponding  to  each  of  the  N^.  spectra 

G,  an  by  dimensional  array  containing  the  geographic  coordinates  for  each  of  the 
X  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

M,  an  by  dimensional  array  containing  the  geomagnetic  coordinates  for  each  of 
the  X  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

v|/,  an  by  «,  dimensional  array  containing  the  nadir  angle  vj/  for  each  of  the  x  n^. 

pixels. 

OUTPUT; 

S  =  [5|,iS'2,...5'jv  j,  the  vector  of  pixel  indices  corresponding  to  each  of  the  spectra. 

Each  element,  S.,  of  the  vector  consists  of  a  pair  of  integers  (m,  n) 
giving  the  column  number  and  row  number  of  the  pixel  contain  the 
foot  of  the  corresponding  field  line. 

A^„,  the  number  of  distinct  subsatellite  pixels.  Generally,  is  less  than  N^. 

P  =  [tJ,  -  ]’  ^  vector  containing  the  distinct  pixel  indices.  P  is  just  S  with 

duplicates  removed. 

Procedure: 

Step  1 :  Locate  boundaries 

Definition:  ^^e  indices  of  the  spectra  that  mark  the 

boundaries  of  the  auroral  precipitation 

Definition: 

Details  of  procedure  TBD. 

Step  2:  Trace  field  lines 
A^.„=0 

For  1=1  to  do 
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Begin  Field  Line  Trace  at  the  position  (>.,.,9,.,  A,) 

[Precise  field  line  tracing  algorithm  depends  on  the  magnetic  field  model  chosen.] 
End  Field  Line  Trace  whenz  =  1 10  km 
Begin  Pixel  Search 

[Precise  pixel  search  algorithm  is  TBD.] 

End  Pixel  Search  when  pixel  containing  foot  of  field  line  is  found 
5,  ={m,n) 
lfi^\  then 

y  =  1,2,.. .1-1  then 

Px.-S, 

End  If 
Else 

N..  =  1 
/]=5, 

End  If 
End  For 

End  Procedure 

END  MODULE  2.  _ _ 


4.3  Sensor  fusion  analysis  of  subsatellite  pixel  data 

Since  the  characteristics  of  the  precipitating  particles  (electrons  and  protons)  are  known 
for  these  pixels,  the  intensities  of  the  four  spectral  regions  (121.6  nm,  125.6  nm,  LBHl,  and 
LBH2)  sampled  by  the  SSUSI  instrument  can  be  predicted.  However,  the  predictions  depend  on 
assumptions  about  the  neutral  composition,  in  particular  the  relative  column  abundances  of  N2, 
O,  and  O2.  Therefore,  any  differences  between  the  predicted  intensities  and  the  observed 
intensities  may  be  interpreted  in  terms  of  departures  of  the  actual  atmospheric  composition  from 
the  assumed  composition. 


MODULE  3:  Sensor  fusion  data  analysis 

PURPOSE:  To  deduce  0/N2  and  02/N2  ratios  by  comparing  the  observed  intensities  with  the 
intensities  expected  from  the  measured  electron  and  ion  spectra. 

INPUT: 

the  number  of  electron  and/or  ion  spectra  taken  during  the  same  time  interval  as 
the  image 

T  =[t,,T2,...Tjv  ],  the  vector  of  times  (UT)  corresponding  to  each  of  the  spectra 
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^  the  vector  of  satellite  geographic  latitudes  corresponding  to  each  of 

the  spectra 

9  =  [9  p92  ’  •  •  •9  /V,  ]» the  vector  of  satellite  geographic  longitudes  corresponding  to  each  of 
the  spectra 

C'"'  =  the  vector  of  electron  spectra  characteristics  as  determined  by 

the  SSJ/5  data  analysis  module.  Each  element  of  the  vector  is 
an  eight  element  vector: 

~  ’  ^Oe.J » Qoej  ’  ^m.j  ’  ^^e.J  •>  ^\e.j » Q\e.j  ’  ^ie.j  ] 

C*/’)  cSv"^],  the  vector  of  proton  spectra  characteristics  as  determined  by 

the  SSJ/5  data  analysis  module.  Each  element  of  the  vector  is  a 
two  element  vector:  Cf  =  [^op.y’^op,/]- 

G,  an  by  dimensional  array  containing  the  geographic  coordinates  for  each  of  the 
"v  ^  ",  pixels*  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

M,  an  by  dimensional  array  containing  the  geomagnetic  coordinates  for  each  of 
the  X  tiy  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

,  an  by  n  ^.  dimensional  array  containing  the  nadir  angle  y  for  each  of  the  n,  x  « 
pixels. 

an  by  dimensional  array  containing  the  121.6  nm  intensity  I^2\6 1®**  of 
the  X  n^,  pixels 

li35f,,  an  by  «,  dimensional  array  containing  the  135.6  nm  intensity  7,35^  for  each  of 
the  X  pixels 

I|BHP  a”  "v  l^y  "v  dimensional  array  containing  the  LBHl  intensity  /lbhi  each  of  the 
X  n^.  pixels 

Ii  bh2'.  an  "v  hy  dimensional  array  containing  the  LBH2  intensity  /lbh2  for  each  of  the 
X  n^.  pixels 

the  number  of  distinct  subsatellite  pixels.  Generally,  7V„  is  less  than  N^. 

P  =  [/^,  Pj.— a  vector  containing  the  distinct  pixel  indices.  P  is  just  S  with 
duplicates  removed. 

S  =  [5'|,5'2,...5jv  ],  the  vector  of  pixel  indices  corresponding  to  each  of  the  spectra. 

Each  element,  S.,  of  the  vector  consists  of  a  pair  of  integers  (m,  n) 
giving  the  column  number  and  row  number  of  the  pixel  contain  the 
foot  of  the  corresponding  field  line. 

OUTPUT: 

N^,  the  number  of  intervals  in  the  normalized  latitudinal  coordinate  jc  defined  below. 
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X  =  ],  the  vector  of x  coordinate  values,  where  jc  is  a  normalized  latitudinal 

coordinate  such  that  at  the  equatorward  auroral  boundary  x  =  =  0 

and  at  the  poleward  auroral  boundary  x  =  =  I 

vector  of  atomic  oxygen  (O)  scale  factors  necessary  to  bring 

the  predicted  intensities  in  line  with  the  observed  intensities. 
Each  element  of  the  vector  consists  of  a  two  element  vector 
corresponding  to  the  two  local  times  at  which  the  satellite  cuts 
the  auroral  oval: 

the  vector  of  atomic  oxygen  (O2)  scale  factors  necessary  to 

bring  the  predicted  intensities  in  line  with  the  observed 
intensities.  Each  element  of  the  vector  consists  of  a  two 
element  vector  corresponding  to  the  two  local  times  at  which 
the  satellite  cuts  the  auroral  oval:  q*®'*  = 

T  =  [TJ,  7^],  a  two  element  vector  containing  the  two  nominal  local  times  at  which  these 

scale  factors  apply. 


Procedure: 


Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 


Np  is  the  number  of  electron  and  ion  spectra  corresponding  to  a  single 
subsatellite  pixel. 

is  an  eight  element  vector  containing  the  sum  of  the  electron 
characteristics  for  the  Np  spectra  corresponding  to  a  single  subsatellite  pixel. 

is  a  two  element  vector  containing  the  sum  of  the  proton  characteristics 
for  the  Np  spectra  corresponding  to  a  single  subsatellite  pixel. 

is  an  eight  element  vector  containing  the  average  electron  characteristics 
for  a  single  subsatellite  pixel. 

is  a  two  element  vector  containing  the  average  proton  characteristics  for 
a  single  subsatellite  pixel. 

I  =  is  a  four  element  vector 

containing  the  observed  intensities  for  a  single  pixel. 

R  =  J  is  a  two  element  vector  containing  scale  factors  for  O  and  O2. 

RADJ  is  a  function  procedure  that  determines  the  best  values  for  R . 

Xs  =  [^.H » 2  ’  •  •  •  V  ]  is  a  vector  containing  the  normalized  x  coordinates  for 

each  of  the  N^  spectra. 

j  is  a  vector  containing  the  O  scale  factors  at  x^. 
is  a  vector  containing  the  O2  scale  factors  at  x,. 


For  i=l,A^„  do 
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sd''*  =  0 

Determine  how  many  electron  and  ion  spectra  correspond  to  this  pixel. 
Fory-l,A^^.  do 
If  5^.  =  i^  then 


k  =  k  +  \ 

sd^’  =  sd^’  +cj.^’ 

End  If 
End  For 


^  “  [^1216 

R=RADj(ci;',C^''',l) 


X;  - 


^pl~‘^eq 

=  Ro 


End  For 

=  INTRP(x,,p'°',x) 
^(0,)  ^INtrp(x^,t|<°^>,x) 


{NOTE:  RADJ  will  probably 
require  additional  arguments} 


End  Procedure 


END  MODULE  3. 


5.  SSUSI  Off-Track  Data  Analysis  Module 


This  module  processes  the  off-track  pixel  data  ingest  by  the  Subsatellite  module. 


5. 1  Data  analysis  of  off-track  pixel  data 

Since  the  characteristics  of  the  precipitating  particles  (electrons  and  protons)  are  known 
for  these  pixels,  the  intensities  of  the  four  spectral  regions  (121.6  nm,  125.6  nm,  LBHl,  and 
LBH2)  sampled  by  the  SSUSI  instrument  can  be  predicted.  However,  the  predictions  depend  on 
assumptions  about  the  neutral  composition,  in  particular  the  relative  column  abundances  of  N2, 
O,  and  O2.  Therefore,  any  differences  between  the  predicted  intensities  and  the  observed 
intensities  may  be  interpreted  in  terms  of  departures  of  the  actual  atmospheric  composition  from 
the  assumed  composition. 

MODULE  1:  Off-track  pixel  data  analysis 

PURPOSE:  To  deduce  electron  and  proton  characteristics  from  the  SSUSI  image  data. 

INPUT: 

G,  an  by  n,.  dimensional  array  containing  the  geographic  coordinates  for  each  of  the 

X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair  )  • 

M,  an  M,  by  «,  dimensional  array  containing  the  geomagnetic  coordinates  for  each  of  the 

X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair  (A,„,,(p„,). 
T,  an  by  m,,  dimensional  array  containing  the  local  time  for  each  of  the  x  pixels, 

an  by  dimensional  array  containing  the  121.6  nm  intensity  7,2,6  for  each  of  the 
X  n^.  pixels 

1,356,  an  by  n^.  dimensional  array  containing  the  135.6  nm  intensity  7,356  for  each  of  the 
X  n,,  pixels 

Ilbhp  an  n,  by  dimensional  array  containing  the  LBHl  intensity  7,  bhi  for  each  of  the 
X  pixels 

Ilbh2  >  an  n,  by  n^.  dimensional  array  containing  the  LBH2  intensity  7lbh2  for  each  of  the 
X  n,  pixels 

7V„  ,  the  number  of  distinct  subsatellite  pixels. 

p(")  ^  a  vector  containing  the  7V„  distinct  pixel  indices. 

N^,  the  number  of  intervals  in  the  normalized  latitudinal  coordinate  x  defined  below. 

X  =  [jCo,x, ,...X5,  ],  the  vector  of  x  coordinate  values,  where x  is  a  normalized  latitudinal 

coordinate  such  that  at  the  equatorward  auroral  boundary  x  =  x^  =  0  and  at 
the  poleward  auroral  boundary  x  =  x^,  =  1 
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the  vector  of  atomic  oxygen  (O)  scale  factors  necessary  to  bring  the 

predicted  intensities  in  line  with  the  observed  intensities.  Each 
element  of  the  vector  consists  of  a  two  element  vector 
corresponding  to  the  two  local  times  at  which  the  satellite  cuts  the 
auroral  oval:  =  [Ti;®’(7;),Tli°*(r2)]. 

^(o,)  ^  the  vector  of  atomic  oxygen  (O2)  scale  factors  necessary  to  bring 

the  predicted  intensities  in  line  with  the  observed  intensities. 
Each  element  of  the  vector  consists  of  a  two  element  vector 
corresponding  to  the  two  local  times  at  which  the  satellite  cuts 
the  auroral  oval :  t]  ^  =  [^ri  ( 7^ ),  t|  ^  ( 72 ) j . 

j(n)  ^  j'j’li)  g  element  vector  containing  the  two  nominal  local  times  at  which  these 
scale  factors  apply. 

OUTPUT: 

ST^,  an  by  dimensional  array  containing  the  electron  energy  spectra  type  for  each  of  the 

^  "v  pixels 

an  by  dimensional  array  containing  the  characteristic  electron  energy  for  each  of 
the  X  pixels 

Qo^,  an  by  n,,  dimensional  array  containing  the  downward  electron  energy  flux  each  of  the 

n^y-n^  pixels 

,  an  by  dimensional  array  containing  the  Maxwellian  characteristic  energy  for 
shifted  Maxwellians  for  each  of  the  x  n^.  pixels 
E(,^„  an  by  dimensional  array  containing  the  characteristic  proton  energy  for  each  of  the 

X  pixels 

Qq^,  an  by  n,,  dimensional  array  containing  the  downward  proton  energy  flux  each  of  the 

n^yn^  pixels 


Procedure: 

Step  1 :  Determine  auroral  boundaries 

Apply  SSUSI  Single  Sensor  boundary  algorithm. 

Apply  boundary  tracing  algorithm  to  produce  ordered  list  of  boundary  pixels. 

If  boundaries  intersect  image  border,  fit  smooth  curve  to  estimate  boundary  location  outside 
of  pixel. 

Produce  vector  and  b^,  to  tabulate  equatorward  and  poleward  boundary  as  a  function  of 
local  time,  . 

Step  2:  Process  pixels. 
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For  i  =  1  to  do 
For  y  =  1  to  Hy  do 

SS.-  =  FALSE 

k  =  0 
Repeat 

k=k  +  l 

If  R  =  Then  SS,  =  TRUE 

If  K  IJ 

Until  55^  =  TRUE  or  A:  =  iV„ 

If  =  TRUE  Then 

This  is  a  subsatellite  pixel,  use  the  measured  electron  and  proton  characteristics. 
Else 

This  is  not  a  subsatellite  pixel,  calculate  the  electron  and  proton  characteristics. 
Ascertain  equatorward  and  poleward  boundaries  for  this  local  time: 

=  INTRP(b^,T„7;,) 

INTRP(b^,T.,7:) 


=INTRP(n'°’,J,is) 

=  INTRP(n'‘’->,x,ls) 

First  Stage  Analysis  (assume  Maxwellian  electron  spectra): 

NONLIN^"  Max  well  ",/2q,/?02  >-^i2i6’A356’^z,s//i»'^tflw2’  ’  0)f 

End  First  Stage  Analysis 
Second  Stage  Analysis: 

Ifa^<0.05a,  Then 
Pure  electron  aurora 

Else  Ifj2o^<20a,Then 

Combined  electron  and  proton  aurora 
If  Qoe  ^  Q.I,re.Ho,d  THCU 

Stage  1  Analysis  OK:  Do  nothing 
Else 

Redo  Stage  1  Analysis  using  monoenergetic  electrons 
NONLIN^"  Mono  ' ,  Rq  ,  Rq^  »  f|216  »  f  I356  •>  ^ LBH\ >  ^ LBHI  ’  Poe’Qoe  ’  ^Op  ’  Qop  ) 

End  If 
Else 

Pure  proton  aurora 
End  If 

End  Second  Stage  Analysis 
End  If 
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End  For 
End  For 

End  Procedure 

END  MODULE  1. 


5.2  Refinements 

In  future  version  of  the  algorithm,  it  may  be  possible  to  refine  the  analysis  by 
interpolating  the  and  functions  using  values  generated  by  the  pure  electron  and  pure 
proton  analysis. 


6.  Auroral  E-layer  electron  density  profile  parameters 

This  module  uses  the  electron  and  proton  characteristics  and  the  neutral  composition  data 
determined  by  the  previous  modules  to  estimate  N^E  and  h^E. 

Ionization  is  produced  by  electron  impact,  proton  impact,  and  photon  impact.  The 
processes  are  linear  so  the  total  production  rate  of  ions  is  just  the  sum  of  the  production  rates 
from  the  three  processes.  The  electron  density  is  obtained  using  a  simple  chemical  model  with 
an  effective  recombination  rate.  Unlike  the  single  sensor  algorithm,  which  used  a  single  altitude 
profde  for  the  recombination  rate,  the  new  algorithm  allows  the  recombination  rate  to  depend  on 
neutral  composition. 

MODULE  1:  Calculation  of  Auroral  E-layer  parameters. 

PURPOSE:  To  calculate  the  E-layer  parameters  N^E  and  h^E. 

INPUT: 

the  number  of  columns  of  pixels  in  the  image. 

,  the  number  of  rows  of  pixels  in  the  image. 

Q^^Jy,\hQ  incident  solar  EUV  flux  (erg  cm'^  s'‘)  as  determined  by  the  SSUSI  single  sensor 
algorithm. 

G,  an  n^  by  dimensional  array  containing  the  geographic  coordinates  for  each  of  the 
n^  X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

M,  an  n^  by  dimensional  array  containing  the  geomagnetic  coordinates  for  each  of  the 
n^  X  n^.  pixels.  Each  element  of  the  array  consists  of  a  latitude-longitude  pair 

T,  an  by  dimensional  array  containing  the  local  time  for  each  of  the  n^  x  n^.  pixels. 

X ,  an  by  «,  dimensional  array  containing  the  solar  zenith  angle  for  each  of  the  n^  x  n^. 

pixels. 
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an  by  dimensional  array  containing  the  atomic  oxygen  (O)  scale  factors  for  each 
of  the  X  n^.  pixels. 

an  by  dimensional  array  containing  the  molecular  oxygen  (Cb)  scale  factors  for 
each  of  the  x  pixels. 

ST^,  an  by  n^,  dimensional  array  containing  the  electron  energy  spectra  type  for  each  of  the 
X  pixels 

an  by  dimensional  array  containing  the  characteristic  electron  energy  for  each  of 
the  «,  X  n^.  pixels 

an  by  dimensional  array  containing  the  downward  electron  energy  flux  each  of  the 
X  n^.  pixels 

,  an  by  dimensional  array  containing  the  Maxwellian  characteristic  energy  for 
shifted  Maxwellians  for  each  of  the  x  pixels 
Egp,  an  by  dimensional  array  containing  the  characteristic  proton  energy  for  each  of  the 

X  pixels 

Qo^,  an  by  dimensional  array  containing  the  downward  proton  energy  flux  for  each  of 
the  X  n^.  pixels 


OUTPUT: 

N^E,  an  by  dimensional  array  containing  peak  E-layer  electron  density  (cm'  )  for  each 
of  the  X  n^.  pixels 

h^E,  an  by  n,,  dimensional  array  containing  the  height  of  the  E-layer  (km)  for  each  of  the 
X  Mj,  pixels 


Procedure: 

Definition: 

Definition: 

Definition; 

Definition: 

Definition: 

Definition: 

Definition: 

Definition: 


Definition: 


is  the  height  (km)  of  the  peak  ionization  due  to  electron  impact. 
hp  is  the  height  (km)  of  the  peak  ionization  due  to  proton  impact. 
hp^  =  1 08  km  is  the  height  of  the  peak  photoionization  rate  at  the  subsolar 
point. 

is  the  peak  ionization  rate  (cm'^  s'')  due  to  electron  impact. 

Pp  is  the  peak  ionization  rate  (cm'^  s  ')  due  to  proton  impact. 

is  a  scale  height  (km)  that  defines  the  shape  of  the  Chapman  function  that 
describes  the  altitude  dependence  of  the  electron  impact  ionization  rate. 

Hp  is  a  scale  height  (km)  that  defines  the  shape  of  the  Chapman  function  that 

describes  the  altitude  dependence  of  the  proton  impact  ionization  rate. 

Hpf,  =  9  km  is  the  neutral  scale  height  near  hp^ .  It  defines  the  shape  of  the 

Chapman  function  that  describes  the  altitude  dependence  of  the 
photoionization  rate. 

N.  is  the  number  of  points  on  the  predefined  altitude  grid  used  in  the 
algorithm. 
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Definition;  z  ==  ,Z2,...z^  j  is  the  predefined  altitude  grid  on  which  the  electron  density 

profile  is  calculated  prior  to  locating  the  peak  density. 

Definition:  is  the  reduced  height  used  in  the  calculation  of  the  electron  impact 

ionization  altitude  profile. 

Definition:  is  the  reduced  height  used  in  the  calculation  of  the  proton  impact 

ionization  altitude  profile. 

Definition;  is  the  reduced  height  used  in  the  calculation  of  the  photoionization 
altitude  profile. 

Definition:  is  the  electron  impact  ionization  rate  at  a  specified  altitude. 

Definition:  is  the  proton  impact  ionization  rate  at  a  specified  altitude. 

Definition:  p^^  is  the  photoionization  rate  at  a  specified  altitude. 

Definition:  n  =  [^«, ,  , . . .  j  is  a  vector  containing  the  electron  densities  corresponding 

to  the  altitudes  z. 

Definition:  a  is  the  effective  recombination  rate  at  a  specified  altitude. 

Definition;  HEIGHT  is  a  function  procedure  that  calculates  the  height  of  the  ionization 
peak  from  the  particle  characteristics.  The  argument  list  depends  on  the  type 
of  particle. 

Definition:  PEAK  is  a  function  procedure  that  calculates  the  peak  ionization  rate  for  an 
energy  flux  of  1  erg  cm'^  s  '  given  the  particle  characteristics.  The  argument 
list  depends  on  the  type  of  particle. 

Definition:  Ch  is  a  function  procedure  that  calculates  the  Chapman  Grazing  Incidence 
function.  The  argument  is  the  solar  zenith  angle.  Ch(0)  =  1 . 

Definition:  CHAPMAN  is  a  function  procedure  that  calculates  a  Chapman  function,  S,  of 
the  form 

5(w,x)  =  exp[l-w-e^"Ch(x)] 

Definition:  RECOMB  is  a  function  procedure  that  calculates  the  effective  recombination 
rate  as  a  function  of  altitude  and  neutral  composition. 

For  /  =  1  to  do 

For  y  =  1  to  do 

=  HEIGHT("  electron", 
h^=  HEI0HT("  proton", 

/;  =  Co,  X  PEAK(”electron",57;, £„„£») 

/>,=a,xPEAK("proton",£„,) 

„  5.25x10'" 


u,,  = 


Zk-K 


H, 


^k-K 

z,  -h 


ph 


H, 


ph 


=  P,  X  CHAPMAN(u„0) 

Pp  =  Pp^  CHAPMAN(w^,0) 

Pph  ~  Qevv  ^  CHAPMAN^Mp;,,%  j 


m  ij 

End  For 
End  For 


jPe+Pp+Pph 

V  a 

End  For 

N,E„  =  MAX(n,*^) 

K.E,:  =  Z. 


End  Procedure 

END  MODULE  1. 
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